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AUTOMATED  DETECTION  TRACKING  SYSTEMS 
FOR  ACTIVE  SONAR 


INTRODUCTION 

The  Information  and  viewgraphs  in  this  document  were  presented  at  a 
wiSC/Code  331  technical  meeting,  a  more  complete  report  will  be  published 
in  the  future. 
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Viewgraph  1 

Automated  detection  and  tracking  systems  for  both  active  and  passive 
sonar  employ  sequential  decision  analysis  to  detect  and  track  potential 
targets.  The  sequential  procedure  also  drops  false  tracks  from  memory. 

This  is  an  important  operation  in  order  to  maintain  finite  computer  memory 
limits. 

The  sequential  procedure  is  a  natural  way  to  detect  a  signal.  Consider 
an  operator  viewing  a  display.  He  will  announce  a  target's  presence  when 
he's  sure  a  target  is  there.  This  is  the  basic  philosophy  behind  sequential 
detection.  A  target  is  detected  only  when  there’s  sufficient  cause.  But 
unlike  an  operator,  potential  tracks  due  to  noise  can  be  discarded  by  the 
sequential  procedure. 

In  this  document  I  will  present  the  fundamentals  of  sequential  detection. 
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Viewgraph  2 

These  are  the  topics  that  will  be  discussed  in  this  document.  The 
sequential  detector  is  compared  with  a  fixed  sample  detector  in  order  to 
show  the  sequential  detector's  optimum  property  of  minimizing  the  average 
detection  time.  The  performance  measure  used  is  the  relative  efficiency.  A 
sequential  detector  with  a  forced  decision  boundary  is  also  discussed. 

Ihese  results  are  important  for  real  tactical  situations. 

The  one-dimensional  and  two-dimensional  Gaussian  cases  are  considered. 
Then  the  Rayleigh  fluctuating  target  model  of  active  sonar  is  discussed  in 
deta  i  1 . 
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Viewgraph  3 

The  sequential  detection  procedure  will  be  applied  at  the  matched  filter 
output.  The  active  transmission  may  be  CW,  LFM,  or  coded  pulses.  The 
likelihood  ratio  is  constructed  at  the  matched  filter  output  to  obtain  the 
optimum  receiver.  Then  the  sequential  analysis  procedure  is  used  to  detect 
potential  target  tracks.  Noise  tracks  are  discarded.  The  advantage  of  the 
sequential  procedure  is  that  it  minimizes  the  average  decision  time.  This 
procedure  may  also  be  important  for  long  ping  durations  in  survei 1  lance .  In 
this  mode  a  target's  presence  may  be  detected  before  the  ping  duration  is 
completed  giving  a  tactical  advantage  to  the  transmitting  ship.  However, 
here  I  wi I  I  confine  my  discussion  to  multiple  transmissions  only.  This  mode 
is  applicable  to  tactical  sonar. 

The  sequential  detection  procedure  was  developed  by  Abraham  Wald. 
Basically,  a  decision  about  the  target's  presence  is  deferred  until  one  of 
two  possible  boundaries  are  reached.  The  log  1  ikel ihood  ratio  output,  Tr, 
for  each  ping  is  compared  with  two  thresholds,  InA  and  1 n B .  If  Tn  >  InA, 
then  a  target  is  present.  If  Tn  <  In8  noise  only  is  present  and  tne 
potential  track  is  dropped  from  memory.  Whereas,  if  1 n 8  <  Tn  <  InA,  a 
decision  is  deferred  until  the  next  ping  is  received.  This  procedure 
minimizes  the  average  detection  time.  However,  the  amount  of  data 
accumulated  is  now  a  random  variable.  In  order  to  prevent  long  decision 
times  from  occasional  ly  occurring  a  non-constant  decision  boundary  is 
sometimes  employed.  This  is  called  a  truncated  or  forced  decision  test, 
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because  at  some  time  a  decision  will  be  forced  to  occur.  Later  I  will  show 
that  if  a  truncated  test  is  used  the  decision  time  can  be  reduced  further, 
but  at  the  expense  of  higher  error  rates.  Nevertheless,  this  may  be  a 
reasonable  trade-off.  The  quantification  of  a  forced  decision  as  given 
later  allows  an  intelligent  trade-off  to  be  made. 

Wald  was  able  to  relate  thresholds  in  terms  of  false  alarm  (a)  and  false 
dismissal  (0)  probabi I ities  by  neglecting  the  excess  over  the  boundaries. 
This  means  that  the  sequential  detection  procedure  terminates  on  the 
boundary.  This  is  a  good  assumption  for  small  s igna 1 -to-noi se  ratio 
problems.  But  for  high  s i gna 1 -to-no i se  ratio  problems  more  exact  methods 
are  needed. 
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Viewgraph  4 

The  performance  of  sequential  detection  procedures  are  usually  based  on 
the  operating  characteristic  function  (OCF)  and  the  average  sample  number 
(ASN)  as  defined  by  Wald. 

The  OCF  is  defined  as  the  probability,  L(h),  of  dismissing  the  false 
track  as  noise  (H^)  as  a  function  of  the  parameter  h.  Whereas,  the 
probability  of  accepting  the  potential  track  as  a  target  (H  )  is,  1  -  L(h). 
This  holds  because  we  assume  that  the  test  will  eventually  terminate. 

The  parameter  h  itself  is  a  function  of  the  signal.  It  can  be  obtained 
from  the  characteri sti c  function  from  the  equation  E[eTnh]  =  1.  This 
equation  will  always  be  satisfied  at  the  points,  h  =  1,  h  =  -1 ,  and  h  -  0. 
However,  the  objective  is  to  find  an  analytic  expression  for  h  which 
satisfies  the  above  equation  for  -1  <  h  <  1. 
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Viewgraph  5 

The  ASN  is  defined  as  the  average  number  of  samples  (pings)  needed  to 
terminate  the  sequential  procedure.  The  equation  for  ASN  (n)  is  shown  here 
for  h*0  and  for  h=0. 


13 


DETECTOR 


RELATIVE  EFFICIENCY 


TD  798b 


Viewgraph  6 

The  next  6  figures  show  asymptotic  results.  Since  it  is  difficult  to 
obtain  general  design  results  for  a  logl ikel ihood  ratio  receiver,  because  it 
depends  on  knowledge  of  the  probability  density  functions,  asymptotic 
performance  results  are  very  useful.  Here  we  assume  that  the  samples 
(pings)  are  statistically  independent  and  that  the  number  of  samples  are 
sufficiently  large  to  assure  that  the  loglikelihood  ratio  approaches  a 
Gaussian  process.  In  this  way  we  can  compare  the  performance  of  a  fixed 
sample  detector  and  a  sequential  detector  using  the  relative  efficiency  as  a 
performance  measure.  However,  non-asymptot i c  results,  which  will  be  given 
later,  correspond  with  these  results  very  well. 

The  relative  efficiency  is  defined  as  the  ratio 

RE  -  N/E(n) 

where,  N  is  the  number  of  samples  (pings)  required  by  a  fixed  sample 
detector  to  achieve  the  desired  a  and  13.  Whereas,  E(n)  is  the  average 
number  of  samples  of  the  sequential  detector  under  the  same  desired  a  and  (3. 
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Viewgraph  7 

The  ASN  is  defined  as  the  average  number  of  pings  needed  to  make  a 
decision  for  prescribed  a  and  (3. 

The  figure  gives  design  curves  as  a  function  of  a  and  (3.  To  generalize 
the  results  the  ASN  has  been  multiplied  by  d,  the  signal-to-noise  ratio. 
Three  conditions  are  plotted  in  the  figure,  h  =  ±  1  and  h  =  0. 
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Viewgraph  8 

In  this  figure  the  relative  efficiency  (RE)  is  plotted  for  h  =  ±  1.  As 
a  and  13  approach  zero  the  RE  approached  6  d8.  This  means  that  the  fixed 
sample  detector  requires  4  times  as  many  pings  as  does  the  sequential 
detector  to  make  a  decision. 

At  a  =  10  ®  and  6  =  10  1  then  the  RE  is  10  dB  at  h  =  1,  and  2  dB  at 
h  =  -1 .  These  results  are  very  significant  and  clearly  show  the  benefit  of 
sequential  detection. 
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Viewgraph  9 

This  figure  shows  the  RE  at  h  =  0  or  the  critical  point.  If  a  =  10 
and  (3  =  101  the  RE  is  still  greater  than  one. 
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Viewgraph  10 

also  important  to  know  what  the  probability  of  detection  (power) 

0.  This  figure  compares  the  power  of  the  sequential  detector  with 
of  the  fixed  sample  detector  at  h  =  0. 
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Viewgraph  11 

In  a  previous  figure  truncating  the  sequential  detector  to  force  a 
decision  at  some  point  in  time  was  discussed.  Here  the  thresholds  are 
modified  by  a  function  G(n).  This  is  a  very  general  approach  to  the  forced 
decision  boundary  problem.  The  performance  measures  are  given  in  this 
figure  for  the  modified  thresholds.  For  the  special  case  shown  truncation 
reduces  the  ASN ,  but  on  the  other  hand  the  error  rates  increase.  The 
derived  mathematical  relationship  can  now  be  used  to  trade-off  the  number  of 
pings  with  the  increased  error  rates.  These  results  also  hold  in  general. 
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Viewgraph  12 

This  figure  shows  a  one-dimensional  Gaussian  example.  Here  the  matched 

filter  output  is  a  unit  variance  Gaussian  process.  The  mean  value  is  e^/2 

under  and  -e^/2  under  Hq.  The  1 og 1 i ke 1 i hood  ratio  is  employed  to 

obtain  T  . 
n 

In  the  next  two  figures  the  results  for  this  example  are  given.  But  the 
data  is  first  quantized  into  two  levels.  This  is  appropriate  since 
information  given  to  operators  on  displays  are  quantized  and  in  general 
computers  require  quantized  data. 
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Viewgraph  13 

_8 

for  the  sequential  detector  were  a  =  10  and 
design  signal-to-noise  (SNR)  was  10  dB. 

The  figure  gives  the  probability  of  terminating  the  sequential  detector 
with  the  acceptance  of  H]  as  a  function  of  n.  The  plot  is  therefore  the 
cumulative  distribution  function  for  n. 

At  10  dB  SNR  Prob(n)  approaches  .9  as  n  increases.  Notice  that  there  is 
a  probability  of  .1  of  terminating  the  test  with  the  acceptance  of  HQ>  but 
this  is  the  desired  error  rate  B. 

If  the  true  SNR  is  less  than  the  design  SNR  then  Prob(n)  <  .9.  This  is 
a  general  result  which  holds  for  all  sequential  detectors. 


The  design  parameters 
(3  =  10  1 ;  whereas  the 
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Viewgraph  14 

In  this  figure  the  design  SNR  was  16  dB.  Now  Prob(n)  approaches  .99 
even  though  a  and  (3  are  the  same  as  before.  The  reason  this  happens  is, 
because,  the  excess  over  the  boundaries  cannot  be  neglected  for  high  SNR. 

As  the  true  SNR  decreases  from  design  conditions  Prob(n)  also  decreases. 
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Viewgraph  15 

Now  the  two-dimensional  Gaussian  example  will  be  considered.  In  this 
case  we  are  assuming  that  the  matched  filter  output  is  correlated  with  the 
previous  ping  output.  This  is  called  a  Markov  process.  The  output  of  a  two 
level  quantizer  will  be  again  considered  in  the  following  data  examples,  but 
now,  instead  of  two  regions  as  in  the  one-dimensional  case,  there  are  4 
quadrants . 
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Viewgraph  16,  17,  and  18 

lhe  next  three  figures  show  the  results  for  the  two-dimensional  Gaussian 
example.  Only  the  correlation  coefficients  were  changed  in  each  figure. 

They  are  .25,  .5,  and  .75,  respectively  for  the  next  three  figures. 

The  previous  figure  and  the  next  three  figures  are  for  the  same  high  SNR 
case  of  16  dB.  However,  as  the  correlation  coefficient  increases  the 
performance  improves  for  the  off  design  SNR  cases.  This  result  may  be 
significant  because  it  represents  a  method  to  improve  active  sonar 
performance.  However,  the  receiver  must  be  designed  from  the  likelihood 
ratio  in  order  to  achieve  improvements.  This  is  true  because  the  likelihood 
ratio  incorporates  correlation  to  achieve  an  optimum  receiver. 

As  far  as  the  figures  are  concerned,  correlation  reduces  the  false 
dismissal  rates  and,  thereby,  increases  detection.  Inspection  of  the 
likelihood  ratio,  on  the  other  hand,  reveals  that  correlation  shifts  its 
weighting  structure  to  favor  detection  when  correlated  data  are  present. 
Whereas,  potential  tracks  produced  by  uncorrelated  noise  data  are  quickly 
rejected . 
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Viewgraph  79 

The  next  example  is  for  the  Rayleigh  fluctuating  target.  This  model  is 
currently  employed  in  active  sonar.  Therefore,  a  complete  formulation  of 
the  sequential  detector  using  Rayleigh  statistics  are  given. 

Here  the  envelope  of  the  matched  filter  output  is  represented  by  x.  The 
probability  density  functions  under  HQ  (noise  only)  and  H  (signal  and 
noise)  are  shown.  The  logl i ke 1 i hood  ratio  from  these  two  densities  is  given 
by  T..  This  gives  the  optimum  receiver  for  the  Rayleigh  fluctuating 
target  model.  From  this  information  and  the  thresholds  the  ASN  can  be 
computed.  Several  cases  are  shown.  For  the  high  SNR  case  (S1  =  16  dB), 
Wald's  formulation  does  not  give  an  accurate  prediction.  This  is  due  to  the 
assumption  of  neglecting  the  excess  over  the  boundaries.  For  lower  SNR 
Wald's  formulation  is  very  accurate.  As  SNR  decreases  more  pings  are 
required  to  detect  the  target  for  the  prescribed  error  rates.  Notice  how 
the  number  of  pings  required  falls  off  as  SNR  decreases.  This  is  a 
nonlinear  relationship.  For  only  a  3  dB  increase  in  SNR  there  is  a 
significant  decrease  in  the  number  of  pings  required  to  detect  the  target. 
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Viewgraph  20 

In  the  previous  figure  the  ASN  for  the  Rayleigh  fluctuating  target  was 
obtained.  Here  the  operating  characteristic  function  is  evaluated  for  the 
Rayleigh  fluctuating  target. 

The  objective  is  to  solve  the  integral  equation  for  the  parameter  h. 
This  integral  can  be  solved.  The  derived  parametric  relationship  for  h, 
however,  is  nonlinear  with  respect  with  S  and  Sy  The  small  signal 
approx imat ion  is  also  given.  This  relationship  agrees  with  the  generalized 
asymptotic  results  given  previously. 

Once  h  is  obtained  the  operating  characteristic  function  can  be 
evaluated  given  the  thresholds.  As  h  varies  from  1  to  -1 ,  L(h)  varies  from 
1  -  a  to  13. 

Some  specific  results  are  shown. 
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Viewgraph  21 

In  this  example  the  data  at  the  matched  filter  is  again  quantized  into 
two  regions.  The  regions  are  separated  by  T  in  the  figure.  When  the  data 
from  a  ping  falls  into  one  of  the  regions  a  value  is  assigned  based  on  the 
probability  of  being  in  this  region  under  and  .  These  values  are 
denoted  by  b^  and  b^  in  the  figure. 

Ihe  next  figure  shows  the  results  of  the  two-level  quantizer. 
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Viewgraph  22 

The  results  of  the  two-level  quantizer  example  is  shown  in  this  figure. 

_ Q  i 

The  design  SNR  is  16  dB  and  a  is  10  and  13  is  10  as  before.  We  see 
that  for  16  dB  SNR  the  sequential  detector  will  achieve  its  design 
conditions.  However,  as  the  SNR  decreases  from  the  design  SNR  the 
performance  falls  off  and  the  desired  false  dismissal  rate  is  not  achieved. 
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Viewgraph  23 

Automated  detection  and  tracking  systems  which  utilize  sequential 
analysis  are  assured  of  minimizing  their  average  decision  Lime.  This  is 
accomplished  by  incorporating  knowledge  of  the  underlying  probability 
density  functions  based  on  the  likelihood  ratio.  The  resulting  automated 
system  acts  much  like  a  human  operator  in  that  the  sequential  detector 
defers  a  decision  until  a  high  level  of  confidence  in  the  target  is 
obtained.  But  unlike  an  operator  target  tracks  which  accumulate  low  levels 
of  confidence  are  dropped  from  computer  memory.  This  is  an  important 
requirement  because  it  allows  additional  potential  targets  to  be  tracked. 

The  automated  system,  however,  is  sensitive  to  parameter  and  probability 
density  function  mismatch.  There  are  several  methods  to  overcome  these 
disadvantages.  These  methods  will  be  discussed  in  future  reports. 


References 

1.  A.  Nut  tall,  "Signal  Processing  in  Reverberation  -  A  Summary  of 
Performance  Capability,"  NUSC  TM  No.  TC-173-72,  30  August  1972. 

2.  R.  Dwyer  and  L.  Kurz,  "Characteri zing  Partition  Detectors  with 
Stationary  and  Quas  i  -Stationary  Markov  Dependent  Data,"  IEEE  Transactions  on 
Information  Theory,  Vol .  IT-32,  No.  4,  July  1986. 


47/48 
Reverse  Blank 


